Breakdown of the Fermi liquid theory in heavy fermion compounds 
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Abstract 

We review the anomalous properties of heavy fermion compounds like CeCu6-iE Au^ or CeMIns close to a zero 
temperature phase transition called a quantum critical point. Anomalous behavior of the resistivity, specific heat 
and magnetic properties observed in various compounds suggests a fundamental breakdown of the Fermi liquid 
theory. Recent measurements on YbRli2Si2, field-tuned through the quantum critical point, provide interesting 
insights on the evolution of the Fermi liquid close to criticality. We discuss the possibility of a kind of spin-charge 
separation at the quantum critical point and make some remarks about its possible link with the underscreened 
Kondo model. 
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1. Experimental overview and theoretical 
insights 

1.1. experiments 

Fermi liquid theory [1], developed by Landau in 
the mid 1950's, provides the foundation for much of 
our current understanding of electron fluids. Accord- 
ing to the Landau Fermi Liquid theory (LFL) , thermo- 
dynamic and transport properties of a metal at very 
low temperatures are described in terms of weakly in- 
teracting fermions, or "Landau quasiparticles" . Lan- 
dau quasiparticles consist of electrons, surrounded by 
a cloud of spin and charge polarization. They share 
the same quantum numbers as free electrons but their 
masses can be strongly renormalized by the back flow 
of the surrounding fluid. A striking example of the ro- 
bustness of the LFL is provided by heavy fermion be- 
havior, where the effective mass of the quasiparticles 
can be hundreds of times greater than that of a bare 
electron. The physical properties of this fluid follow cer- 
tain characteristic power-law dependences on temper- 



ature. For example, the specific heat has a linear tem- 
perature dependence, the resistivity a quadratic one. 

A wide body of experimental results show that when 
heavy fermion materials are tuned to a zero tempera- 
ture antiferromagnetic transition, their physical prop- 
erties depart qualitatively from LFL theory. There is 
a growing list of heavy electron materials that can be 
tuned into the quantum critical point by, either alloy- 
ing, such as CeCuQ- x Au x [2] or the layered materials 
CeMIns[3], through the direct application of pressure, 
as in case of Ce/ri3[4] or CePd2Si2[5], or via the appli- 
cation of a magnetic field, as in the case of Y&i?/i2<Si2 [6]. 
At the quantum critical point (QCP) the key proper- 
ties suggesting a breakdown of the Fermi liquid are: 

- the specific heat coefficients diverge at the 
QCP[7,8,9]. For most cases the divergence displays 
a logarithmic temperature dependence 
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This implies that the the quasi-particle mass di- 
verges and the Fermi temperature goes to zero at a 
QCP. 
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a quasi-linear temperature dependence of the resis- 
tivity, contrasting with the quadratic dependence of 
the FL [5,6,10]. 
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with e in the range of — 0.6. Several compounds 
including YbRli2Si2 [11] and CeCue-^Aua; [9] have 
perfectly linear resistivity, a property reminiscent of 
the normal phase of high temperature superconduc- 
tors. 

the spin susceptibilities acquire anomalous expo- 
nents 
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with a < 1 for CeCus.gAuo.i, YbRl^Sii-xGea,^ 
(x — 0.05) and CeNiaGe2[5]. In the case of 
CeCu6-iAu2, [9] magnetization measurements re- 
veal that at finite fields H , the differential magnetic 
susceptibility exhibits exhibits H/T scaling. Neu- 
tron scattering measurements on the same material 
[9] reveal the presence of ui /T scaling in the dynamic 
spin sysceptibility, which may be cast into the form 



X _1 (q,w) =/(q) + (^ + T) a . 
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where the function /(q) vanishes in the vicinity of 
the ordering wave vector. The q— independence of 
second term suggests a local origin to the damping of 
the critical spin fluctuations, an observation that has 
stimulated recent efforts to develop a self-consistent 
locally quantum-critical model of the heavy electron 
QCP. [12] Remarkably, a single exponent a ~ 0.75 
governs both the H/T and u>/T scaling. The pres- 
ence of the anomalous exponent a suggests that the 
interaction amongst the quantum critical modes has 
renormalized to strong coupling. This fact is made 
clear by the scaling laws in u)/T and H/T. This type 
of "naive" scaling property (quantities scale accord- 
ing to their dimensions) is the hallmark of a system 
lying below its upper critical dimension. [13] 

1.2. Key theoretical insights 

The experimental observations at criticality sug- 
gest an underlying universality in the physics. [13]. 
The ability to express the dynamical susceptibility of 
Cus.gAuo.i in terms of a scaling function F (w/T, B/T) 
and the observation that the effective Fermi tempera- 
ture T£.renormalizes to zero at the QCP, both suggest 
that underlying physics at the QCP is scale invariant. 
In this situation, we expect that the physical proper- 
ties will be independent of the microscopic scales of the 
material- they are universal. 



Moreover, the observation of anomalous exponents 
in the magnetic susceptibility, the transport and ther- 
modynamic properties are evidence that the quantum 
critical physics may be described by an effective the- 
ory whose coupling constants flow to strong coupling. 
In particular the effective critical lagrangian must lie 
below its upper critical dimension. 

1.3. Spin density wave scenarios 

The constraints at criticality enable us to rule out 
the 3-D spin density wave (SDW) theoretical sce- 
nario [14,15] for the heavy electron quantum critical 
point (figure l,left). In this scenario, magnetism devel- 
ops by the spin polarization of the Fermi surface and 
non-Fermi liquid (NFL) behavior results from Bragg 
diffraction of electrons off quantum critical spin fluctu- 
ations in the magnetization. As shown in figure 1, three 
dimensional spin fluctuations only couple strongly to 
electrons along "hot lines" of the Fermi surface, sep- 
arated by the wave vector Q of the antiferromagnetic 
order. The term "hot lines" refers to the k points on the 
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Fig. 1. Comparison between the 3D SDW Fermi surface topol- 
ogy (left), where hot lines are separated by the scattering wave 
vector and the 2D SDW scenario where a finite part of the 
Fermi surface is "hot" at T — 0. 

Fermi surface which are coupled the magnetic fluctua- 
tions at zero temperature. The width of the hot lines in- 
creases with temperature like yff/T). On the hot lines, 
the electrons are strongly perturbed by the magnetic 
fluctuations and their life time is strongly suppressed. 
More precisely in three dimensions the imaginary part 
of the self-energy varies like Eh ~ w. Their effective 
mass experiences a logarithmical increase with temper- 
ature. Away from the hot lines however, the scattering 
of the electrons off the critical modes is strongly sup- 
pressed - the magnetic fluctuations getting massive. 
Within the 3D SDW scenario, the hot lines are not 
wide enough to destabilize the FL theory. In this situ- 
ation, the critical contributions to the free energy and 
the specific heat coefficient are given by 



F ~ T 5/2 and C/T = 70 - VT 



(6) 



When the temperature is driven to zero at criticality, 
the system eventually flows towards the FL fixed point, 
keeping the universal exponents of the FL theory. 

In order to account for the properties at critical- 
ity, one can imagine that the spin fluid is magnetically 
frustrated, dividing into layers of decoupled plans. In 
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this scenario, the 2D SDW [16,17], the dimension of 
the spin fluid (2D) is lower than the one of the con- 
duction electrons (3D). Within this assumption a fi- 
nite region of the Fermi surface is hot, as shown in 
figure 1, right. This intuitively corresponds to an "in- 
finite" width of the hot lines. This situation- a fully 
anisotropic magnet-, while difficult to realize experi- 
mentally, is strong enough to destabilize the FL the- 
ory, accounting for the logarithmic divergence of the 
specific heat coefficient with temperature as well as 
the linear in T resistivity. Two important remarks are 
due here. First, this scenario is remarkably sensitive 
to magnetic anisotropy. If the 2D spin fluid retains a 
small inter-plane coupling, then Fermi liquid behavior 
will recover at temperature scales smaller than the in- 
terlayer coupling. The second comment relates to the 
upper critical dimension. Close to an antiferromagnetic 
instability at zero temperature, the system experiences 
strong fluctuations in time which scale with respect to 
the length like r ~ £ z . z is the dynamical exponent 
of the QCP, and is equal to 2 for the antiferromagnet. 
The effective dimension of the system is now D = d+z, 
where d is the number of spacial dimensions for the 
spin fluid. In the 2D SDW scenario, D = 2 + 2 = 4so 
the effective lagrangian lies at its upper critical dimen- 
sion, implying marginal coupling constants. This fea- 
ture wouldn't be enough to account for the anomalous 
exponents observed in the magnetic susceptibility. 



2. A study of YbRh 2 Si2: the breakdown of the 
2D SDW scenario 

In the next section, we would like to discuss new 
insights into the the mechanism by which the heavy 
quasiparticle may disintegrate at the QCP, revealed 
by recent studies of the compound YbRh2(Sii-a; Ge^ 
[18,19]. The undoped (x = 0) compound lies remark- 
ably close to a QCP, with a tiny antiferromagnetic or- 
dering temperature Tjv = IQmK [11] that can be sup- 
pressed to zero by a small magnetic field of just B c — 
0.06T perpendicular to the easy axis. By allowing tiny 
amount of Germanium into the crystal, it is possible 
to fine-tune the Neel temperature Tjv and the critical 
field B c to zero, without appreciably increasing the dis- 
order in the material. The materials can then be con- 
trollably driven back into the FL state by application 
of a small magnetic field. 



coefficient 7 shows a strong upturn below 30mK (figure 
2). The electronic nature the upturn has been checked 
first by carefully removing the nuclear contribution, 
second by the observation that under a small magnetic 
field, the electronic contribution in the FL phase grad- 
ually saturates the upturns, as one gets closer to the 
QCP. Such an upturn of electronic nature reveals a 
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Fig. 2. Strong upturn in the specific heat coefficient of 
YbRh 2 (Sii_ x Gc x ) 2 for x = 0.05 [19]. One observes how the 
upturn fills in as one reaches the QCP by decreasing the ap- 
plied magnetic field. 

strong departure from any SDW scenario, where the 
maximum possible increase of the specific heat coeffi- 
cient is logarithmic. More specifically, the behavior of 
7 can be cast into the form 

j(T,B)=B- a ^) (7) 

with a ~ 0.33 and r) ~ 1. $ is a scaling function of 
the temperature and the magnetic field. The form of 
(7) is remarkable in many respects. It reveals that two 
energy scales govern the behavior of the electrons in the 
FL phase. A first energy scale 2> is associated to the 
inverse of the density of states. The second energy scale 
To(B), characteristic of the excitations of the system 
appears in the scaling function 3>. These two energies 
vary with the applied magnetic field as follows: 

Tf ~ B a and Tq(B) ~ B (8) 

The scaling form (7) of 7 ensures that at zero magnetic 
field 7(T,S = 0) oc T~ a/,? . Since the exponent a is 
smaller than one, no residual entropy is sitting at the 
QCP. 

2.2. Kadowaki- Woods ratio is weakly field dependent 



2.1. The upturn in the specific heat 

Several unusual features are revealed at criticality. 
In addition to the previously observed logarithmic in- 
crease when lowering the temperature, the specific heat 



An alternative way of investigate the scattering be- 
havior at the QCP is to approach it from the FL phase, 
gradually reducing the applied magnetic field towards 
its critical value B c . Away from criticality, the resis- 
tivity is quadratic Ap(T) = p(T) - p = A(B)T 2 (p 
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is the residual resistivity), below a characteristic tem- 
perature Tp. The temperature at which the T 2 behav- 
ior, characteristic of the Landau FL theory appears, 
sets the crossover between the heavy fermion regime 
and the quantum critical regime. Inside the FL phase, 
even though the quantum fluctuations are not strong 
enough to destabilize the FL, the electron undergoes 
virtual scattering off the critical magnetic modes. The 
variation of A(B) close to the QCP provides insight 
into these virtual excitations. Initial measurements on 
germinium doped YbRh^Sii-xGex show that A{B) ~ 
1/(5 — B c ) in the approach to the QCP. Over a wide 
range of field, the Kadowaki Woods ratio is found to 
be approximately constant, with a value 



K = A/ {'to) ~ 5-8fincmK 2 mof/J 
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is found to agree within 40% with the empirical Kad- 
owaki Woods ratio. More recent measurements closer 
to the QCP show a slow upturn in K at very small val- 
ues of B — B c . [19] Constancy of the Kadowaki woods 
ratio implies that the transport scattering rate scales 
strictly with the square of the renomalized T F (B) of 
the heavy electron fluid. A truly field independent Kad- 
owaki Woods ratio would indicate that the momentum 
dependence of the scattering amplitude does not renor- 
malize with T F (B). 



specific heat are not involved in the charge transport. 
What is the origin of this separation? 

One of the challenges here is to understand the 
clear separation between the single ion Kondo temper- 
ature Tk , which clearly remains finite at the QCP, and 
the renormalized T F (B) which is driven continuously 
to zero. Burdin et al.[20] have suggested the analogy 
with pre-formed pairs in a superconductor, according 
to which Tk is the temperature at which pre-formed 
local Kondo singlets develop and T F is the scale at 
which phase coherence sets in to form mobile quasi- 
particles. By analogy one can consider the QCP in 
heavy fermions as a king of preformed scenario for the 
Kondo effect. In such a picture, the presence of very 
strong antiferromagnetic fluctuations close to a QCP 
dephases the heavy fermion composite quasi-particles 
at increasingly low temperatures, as one gets closer to 
the QCP. Eventually at the QCP the coherence tem- 
perature is driven to zero, leading to the disintegration 
of the heavy composite electron. 
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2.3. The 2D SDW scenario is ruled out 

This set of data allows us to rule out the 2D SDW 
scenario. While the 2d SDW scenario leads to the re- 
sult A oc 1/{B — B c ) as observed experimentally, the 
soft 2D spin fluctuations produce only a small renor- 
malization in the heavy electron density of states with 
7 oc In 1/{B — B c ). This would give a strongly diver- 
gent Kadowaki Woods ratio in the approach to the 
QCP: Ksdw oc 1/ ({B - B c )]n 2 {B - B )). The weak 
field dependence of the Kadowaki Woods ratio suggests 
that the most singular quasiparticle scattering ampli- 
tudes have a far weaker momentum dependence than 
expected in a spin density wave scenario. 



3. The breakup of the heavy electron- is there 
a preformed Kondo effect at the QCP? 

The study of YbRh2(Si2- :c Ge a: ) reveals a very un- 
usual property. While the specific heat coefficient ex- 
periences a strong upturn below T — 30mK, there is 
no reflection of this sudden loss of entropy in the trans- 
port properties: the resistivity stays linear in T for 
the whole temperature regime. This suggests that the 
degrees of freedom associated with the upturn in the 
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